BLOW-UP IN HIGHER-ORDER REACT ION-DIFFUSION 
AND WAVE EQUATIONS: HOW ^loglog FACTOR OCCURS 



V.A. GALAKTIONOV 

Abstract. The origin of non self-similar blow-up in higher-order reaction-diffusion 
(parabolic) or wave (hyperbolic) equations with typical models 

u t = u 2 (-u xxxx + u) or u a — u 2 (-u xxxx + u) in (-L, L) x (0, T), 

with zero Dirichlet boundary conditions at x = ±L, where L > Lq G (§i 7r ) is discussed. 
The rate of blow-up is shown to get an extra universal ^/ln | ln(T — t)\ factor in addition 
to the standard similarity one ■ The explanation is based on matching with the 

so-called logarithmic travelling waves as group invariant solutions of the equation. 

Some links and similarities with double-log blow-up terms occurring in earlier studies 
of plasma physics parabolic equations and the nonlinear critical Schrodinger equation 
are discussed. On the other hand, y/\n | ln(T — t)\ obtained in Petrovskii's boundary 
regularity study for the heat equation in 1934 was the first its appearance in PDE theory. 



1. Introduction: THREE mysterious blow-up v^og log of PDE theory 

1.1. On our main goal. In the middle of the 1980s, in the study of singularity formation 
phenomena of blow-up in a reaction-diffusion equation from plasma physics and, almost 
simultaneously, in self-focusing for the cubic nonlinear Schrodinger equation, the physical 
and formal asymptotic methods, to say nothing about rigorous justifications, faced an 
extremely difficult issue of appearance the so-called double logarithmic, or i/log log factor: 

(1.1) Ao(t) ~ y/]n\ki(T -t)\ as t^T~. 

Here T < oo is the blow-up time in the sense that the solution u = u(x, t) of the PDE 
under consideration is well-defined and is classic for all t e (0, T), but gets unbounded^: 

(1.2) lim^T- sup x \u(x, t)\ = +oo. 

The goal of this paper is to introduce a number of higher-order scaling invariant nonlin- 
ear PDEs of parabolic, hyperbolic, and nonlinear dispersion types, which can exhibit the 
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factor (11 .ip in their blow-up behaviour. In fact, we are going to address a wider question: 



1.3) under which assumptions on PDEs, (II. ip occurs in blow-up or not, 



and how the absence of fll.l[) affects the generic blow-up behaviour. The mysterious 
(becoming not that much after a proof is found) blow-up factor (II. ip is already well- 
known for a few PDEs, so we are inevitably obliged to begin with this amazing history of 
the twentieth century. 

1.2. FIRST Vlog log in classic parabolic theory: boundary regularity and Petro- 
vskii's criterion, 1934. It is truly amazing that its actual origin lies in the heart of PDE 
theory: regularity of a boundary point. For the Dirichlet problem for the Laplace equa- 
tion, this study, began by Green, Gauss, Lord Kelvin, Dirichlet in the first half of the 
nineteenth century and by many other great mathematicians, was completed by Wiener 
in 1924 [62], who derived his famous regularity criterion (a necessary and sufficient con- 
dition). A detailed history of potential theory can be found in Kellogg [311 pp. 277-285]. 

The same regularity question for the heat equation in a non-cylindrical domain was in 
1934 initiated by Petrovskii [471148] . where the double-log actually occurred for first time. 
This is the question on irregular or regular point (x,t) = (0,T) for the ID heat equation 



Ut = u X x in Qt = {\ x \ < R(t), < t < T}, where R(t) — ► + as t — ► T 
with bounded smooth data u{x, 0) = u (x) on [— i?(0), -R(O)]. 



Here the lateral boundary {x = ±R(t)} is given by a function R(t) that is assumed to be 
positive and C^-smooth for all < t < T and is allowed to have a singularity of R'(t) at 
t = T~ only. Then the value of u(x, t) is studied at the end "blow-up" point (0, T), to 
which the domain Qt "shrinks" as t —>■ T~ . Thus, (x,t) = (0,T) is regular, if any value 
of the solution u(x, t) can be prescribed there by continuity as a standard boundary value 
on 8Qt- Otherwise the point is irregular, if the value u(0,T) is not arbitrary and is given 
by the evolution as t — > T~. 

Petrovskii in 1934-35 [47] 148], in particular, established the following: 



1.5) 



(i) R(t) = 2^/T^t y/\n\]n(T-t)\ (0, T) is regular, and 

(ii) R{t) = 2(1 + e)VT^t </\n | ln(T - t)\, e>0 (0, T) is irregular; 



see the most recent survey in [25] for a full list of references. As far as we know, (11.51) this 
is the first clear appearance of the blow-up y/\og log factor in PDE theory. 



1.3. SECOND Vlog log in the NLSE (1985): the origin and beginning of the 
log-log story in nonlinear PDE theory. We now return to the true blow-up problems, 
but will begin with another famous PDE, which is not a subject of the present paper, 
but used to be well-recognized as the source of the log-log. It has been well-accepted that 
the origin (besides Petrovskii's result more belonging to regularity and even probability 
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theory) of such an unusual rate of blow-up including the ->/ln | ln(T — t)\ factor is the 
nonlinear focusing Schrddinger equation (the NLSE) with the critical exponent 

(1.6) iut = -Au - \u\ p - x u in R N x (0,T), where p = p = 1 + ±, 

which is a fundamental model of water wave theory (N = 1), nonlinear optics (mainly 
N = 2), and plasma physics (TV = 3, as a limiting case of Zakharov's model of Langmuir 
waves, 1972). It seems that G.M. Fraiman in 1985 [20] (see also [55], and a full history of 
asymptotics in [571 P- H5]) was the first, who formally derived the log-log correction to 
the blow-up rate for the cubic case in dimensional N = 2. An alternative way of derivation 
of the log-log was presented in [321 H3] (1988). See the monograph [57] for a full history 
and extra details, and [IB] for an extended list of further qualitative results. These first 
results were based on formal ideas saying that the L°°-blow-up rate for some classes of 
solutions with L 2 -norm slightly above that for the ground state Q(x), satisfies 

N 

(1.9) sup x \u(x, t)\ ~ 

We refer to a number of more recent papers, where first qualitative and formal estimates 
obtained a proper mathematical justification |3 [TDl IIBl HSl EH ESI ESI SDl SD SSI S3 E2] . 
Many questions remain still open in view of the general complexity of the NLSE (II. 6p . 
especially in the mult i- dimensional geometry, with N > 2, where even the representation 
(II. 9p is questionable (other norms can be more appropriate). 

In general, blow-up such as (11.91) for the NLSE (11.61) with the strong L 2 -conservation 
property (the dynamical system being symplectic-Hamiltonian), can be characterized as 

2 There are two simple misprints in |20j (in both Russian and English versions that are identical) that 
are convenient to know: (i) the formula (5.8) must read 

(1.7) "« - (1 - |C | 2 )P cr /a 3 £ c V (instead of «fi = (1 - \C \ 2 P CI ) / a 3 , ")• 

In the process of its derivation, three lines above (5.8) the reference "(3.25)" (a nonexistent formula) 
seems should be replaced by "(3.17)". Then (5.14) implies (5.15), i.e., "| C 1 2 — 1 ~ l/ln 2 r", where 
r ~ — ln(to — i), so (|1.7p yields the result: 

( L8 ) in Wffl => ffl (*)~ VET^Sji" aS t->t °' 

with, according to (3.1), being the desired extra non-self-similar blow-up factor as in (|1.9j) for N = 2. 

(ii) Thus, according to the correct asymptotics in (|1.8[) . in the final formula (5.16) in [20 (cf. p. 400 
in the Russian version), the second In was incidentally missing, as earlier (and for the first time?) noted 
in [57l p. 116]. Overall, Fraiman's derivation of log-log looks very solid and formally well-justified from 
the point of view of standard asymptotic methods; in particular, it does not use a "homotopying" in the 
dimension N — > 2~ , which was used in some other papers a few years later, until perfect solution via 
solid operator-functional analysis approaches developed by Perclman, Merle, Raphael, and others (see 
references below) in the twenty first century only. 

A single log-correction was suggested earlier by V.I. Talanov in the 1970s (see [59]) and by D. Wood in 
the 1980s (see also references in [35.); V.E. Zakharov also claimed earlier derivation of the double log-log 
[63] . see some details and results of the Russian School in [33]. Self- focusing (blow-up) itself in the NLSE 
was proclaimed in the middle of the 1960s, [TTI |3"U1 12], 
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ln|ln(T-t)| 
T-t 



as t 



corresponding to a centre- sub space behaviour relative to the ID manifold of the rescaled 
ground states {Qoj(x)}, though sometimes this is again not well- understood and question- 
able (but this is of importance for other models with a different reason for the factor). 

1.4. THIRD \/log log: a cubic second-order reaction-diffusion model (1985). 

As rather customary in PDE theory, the above log-log factor was almost simultaneously 
(with Fraiman's research [2U]) detected in reaction-diffusion theory, and was involved 
into mathematical blow-up area by Friedman and McLeod in 1986 [22], who studied the 
following parabolic equation from plasma physics: 

(1.10) u t = u 2 {u xx + u) in (—L,L) x (0,T), u(±L,t) = 0, where L > \ . 

More precisely, according to B.C. Low [37_J1, f ll.lOp can be considered as a model of resistive 
diffusion of a force-free magnetic field in a plasma. The multi-dimensional version of the 
model assumes replacing u xx by the Laplacian Au. Actually, Low proposed the original 
nonstationary model, in which a contraction of the magnetic field becomes infinite after 
a finite time provided that the free-force condition is always satisfied. 

It is worth mentioning here that plasma physics with magnetic fields is a permanent 
source of various reaction-diffusion problems, where evolution of "frozen" magnetic field 
is typically described by fast-diffusion-like and other quadratic ambipolar mechanisms. 
For convenience and better understanding a hierarchy of the corresponding parabolic 
models, we present a typical multi-term RD-type equation for the magnetic field B [9J: 
B( = Vx (v n x B + vd x B — AV x B), where A is the ordinary magnetic diffusivity, 
and vd = vj — v n is the drift velocity between ions and neutrals, which is proportional to 
the Lorentz force, vd ~ ( -^^ xB (this is the mechanism of cubic nonlinearities). Overall, 
this yields the following equation: 

B t = Vx[v n xB + ™B-(A + A AD )VxB], where A AD = 

is the ambipolar diffusion coefficient. In addition, there can appear the related ID complex 
cubic PDEs of the PME-type (i.e., unlike f ll.lOp . with a divergent diffusion operator) 
Bt = (^(\B\ 2 )z) z (without ordinary magnetic diffusion, no neutral velocity, and very 
small ionization fractions), which creates sharp structures in finite time, [9, p. L92]. 

First qualitative results for fll.lOp including conjecture of blow-up were obtained in 
[60~t |6lJ^. On the other hand, (jl.lOp appears in curve-shortening flows by curvature[3]. 
Namely, the equation X t = kN of the evolution of a curve X : S 1 x [0, T) — ► M 2 on the 
plane, where k is the curvature and N is the unit normal, is reduced to 

k t = k 2 (k ee + k) on S 1 x (0, T) (S 1 = {0 < 9 < 2ty} is the unit circle in R 2 ) . 



The plasma and current sheet collapse is in fact forming a singularity, suggested by B.C. Low...", 
[551 p. 3]. 

4 The collapse "..has been predicted analytically by Watterson (1986), who used a simplified one- 
dimensional, force-free model to show that a radially decreasing magnetic field profile with a reversal of 
the axial field cannot be maintained indefinitely", [S51 p. 508]. 
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On derivation, history, blow-up patterns, and other details, see papers [31 0] to be referred 
to again. For earlier results and blow-up similarity solutions for such a flow, see e.g., pQ. 

First, Friedman and McLeod [22, § 2] proved blow-up of solutions if L > | by using 
standard at that time barrier and Maximum Principle techniques. Among more delicate 
results, they proved a rather rare (again at that time) result on the blow-up set of the 
solution u(x, t): 

(1.11) B[uq] = {iG (— L, L) : 3 {x n } — > x, {t n } — > T~ such that u(x n ,t n ) — > +00} . 
By the definition, B[uq] is closed. Then, for smooth symmetric monotone data [221 § 4] 

(1.12) u (-x) = u {x), u' (x) < for x > B[u ) = [-§, |]. 

The next question is then the description of the blow-up behaviour in B[uo] as t — > T~, 
which turns out to be a difficult issue. Namely, the authors showed by some involved 
mathematics associated with the equivalent integral equation, obtained by dividing by u 2 
and inverting the operator D 2 + I on the right-hand side, that the blow-up behaviour 
does not obey the dimensional similarity law associated with the separation of variables: 

(1.13) u(x,t) ~ t=0(z) as t^T => \e = 6 2 {6" + 6). 

More precisely, unlike (I1.13p . they showed that the behaviour is governed by the non 
self-similar law 

(1.14) u(x,t) ~ cosx, where A(t) — > +00 as t — > T _ . 

In addition, it was mentioned on the last page [221 p. 80] that in Watterson's Thesis in 
198^3 [60] the ([TID ansatz was madcl- An involved proof of (11 .11) and the behaviour (11.141) 
was achieved later on in 0] with a high sharpness: the result therein was 

(1.15) A(t) = vTn I ln(T - t)\(l + o(l)) as t -»• T~ 

(we will explain why the multiplier 1 • ^77. always occurs all the time here on matching). 

It is convenient for further use to discuss the origin of the structure (11.141) . The cosx 
is obviously a formal stationary solution of (11.101) . 

(1.16) f(x) = cosx: fxx + f = 0, 

which nevertheless does not satisfy the Dirichlet boundary conditions in fll.lOf) . since 
L > |. Therefore, (11. 14j) can be interpreted as an evolution as t —>■ T~ along the ID 
manifold {/icosx, \i > 0} of stationary solutions, where the resulting slow-growing factor 
(II. ip represents a kind of centre manifold (or subspace) behaviour corresponding to the 
linearized rescaled operator. Of course, the crucial argument establishing (II. II) is fully 
dependent on the matching of the blow-up structure (11.141) with bounded solutions near 
x = ±L through the special behaviour at the localization end points x = ±|. 

5 I.e., in the same 1985 year as Fraiman's conclusion for the cubic NLSE [20]! What nonlinear theory 
can explain this historical "time-evolution" coincidence? 

6 As explained, "... on the basis of numerical evidence", which is hardly believed even in the twenty- 
first century: log-log in blow-up is extremely difficult to catch numerically since a huge solution growth 
of hundreds of orders is necessary. 
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1.5. On a "topological" matching: main idea. For further applications to more diffi- 
cult nonlinear PDEs, we need to explain the qualitative origin of our rather "topological" 
asymptotic analysis. We claim that the actual origin of the log-log factor ( 11. ip in the 
RD model (11.101) can be associated with the existence of the log-travelling wave solutions 
(log-TW) of the equation (11.101) of the following form 

(1.17) u log {x,t) = ^= t 9(v), r ? = x + Aln(T-t) (A ^ 0). 

Substituting (071) into (QUI) yields the following ODE: 

(1-18) \g-\g' = g\g" + g). 

These solutions can be obtained as the result of the invariance of ( 11.101) relative to a one- 
parameter group of scaling-translating transformations; this was first proved by Ovsian- 
nikov in 1959 [S]. 

We show that similarity solutions (I1.17P play a key role at the end points x = ±| of the 
localization domain x G (— f , §), and serve as a "transitional mechanism" from the Inner 
Region-II of bounded solutions for |x| > f into the internal Inner Region-I {\x\ < |} with 
the blow-up behaviour (jl,14p . This idea is explained in Figured], which will be used in 
greater detail below. Since in some asymptotic ^-interval, solutions g(rj) of (11.181) have 
an extra logarithmic correction in the asymptotics, i.e., 

(1.19) g{rj) ~ (-r})y/ln(-r}) for r] < -1 

(in this intermediate region, the term g 3 in (11. 18ft ought to be neglected), the combination 
of two logs: 

(i) y/\n(-ri) in f[TTT9|) . and 

(ii) ln(T — t) in the //-variable in (I1.17p . 

lead, on a special formal matching, to the double logs in (ll.ip in the RD model (ll.lOlf l. The 
later difficult proof such a blow behaviour with the factor ( 11. ip in Angenent-Velazquez 
jl] of the behaviour (I1.14p . (II .ip essentially uses the rescaled variables corresponding to 
the log-TWs moving frame (11.171) and other related issues. 

In what follows, to avoid very complicated technical calculus and further justification 
(which are actually nonexistent for higher-order PDEs under consideration), we prefer to 
keep the ideology of the above "topological" matching with a clear geometrical meaning 
to be explained in Figure [1] in Section [2j 

1.6. No log-log for the divergent RD model with finite propagation. The corre- 
sponding cubic RD model with the divergent diffusion operator has the form 

(1.20) u t = {u 3 ) xx + u 3 in Rx(0,T). 

Answering the question (11.31) . the generic blow-up of nonnegative solutions is then de- 
scribed by the following explicit Zmitrenko-Kurdyumov solution [54"] : 

(1.21) u s (x,t) = -jL=6(x), where \ 6 = (9 3 )" + 9 3 in R, so 

"This idea was developed by the author in discussions with Herrero and Velazquez in Dpto. de 
Matematica Aplicada, Unversidad Complutense de Madrid, 1992, as explains in a few lines in [531 P- 308]. 
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1.22) 9( X ) = l^ COS ^ f ° r \ X \^f 



for \x\ > ^f. 

Concerning the asymptotic stability and other generic evolution properties of the par- 
ticular solution fll.2ip . see [521 Ch. 4] and references therein. In other words, the weak 
solution fll.2ip . fll.22p shows that the blow-up in localized in the domain (— ^f, which 
the heat does not penetrate from for all t 6 (0, T). The rate of blow-up in ( 11.211) is purely 
self-similar and does not contain any extra log- log or other factors as in (11 .11) . 

Indeed, returning again to (11.31) . this absence of the log- log factor can be directly 
connected with the solvability of the ODE in (ll.2ip . which admits a good weak compactly 
supported profile 9(x) given in (11.221) . This profile is suitable for both the Cauchy problem 
in R x (0, T) and for the IBVP in (-L, L) x (0, T) provided that 

(1.23) £>T- 

If (11.231) is violated, other boundary conditions at the end points x = ±^r can be imagined 
that can lead to extra blow-up factors. Though such posed IBVP for (11.201) are typically 
looking rather artificial unlike the above and given below more natural parabolic models. 



1.7. Main models and results: explaining log-log in higher-order RD models 
via transition by log-TWs. We claim that the same blow-up factor as in (II .ip exhibits 
more universality and occurs in other nonlinear scaling-invariant cubic parabolic models 
of higher-order. Moreover, we also claim that, similar to (11.171) . 

y/log log is generated by the log-TWs transition mechanism. 

As a key example, we use the fourth-order RD equation, as a natural extension of 
(11.101) . with zero Dirichlet boundary conditions: 

q 241 (ut = u 2 (-u xxxx + u) in (-L, L) x (0, T), 

\ u — u x — at x = ±L (L > L ). 

Of course, the problem for this higher-order parabolic equation loses any traces of order- 
preserving, comparison, the Maximum Principle, and barrier features that are always key 
ingredients in the mathematical study of the second-order parabolic equations. Moreover, 
even existence-uniqueness results (local) for the degenerate equations such as (11.241) are 
not properly settled. Nevertheless, we put a blind eye on those difficulties, and concen- 
trate on predicting the log-log factor in the regional blow-up behaviour. Actually, there is 
some space to avoid such local difficulties: for instance, we can consider uniformly strictly 
positive solutions with such data on the boundary that the similarity blow-up such as 
(11.211) is impossible (but then some extra speculations are necessary). For positive solu- 
tions u(x, t) > 5o > 0, classic parabolic theory [T|J[2T] applies to guarantee local existence 
and uniqueness of smooth (moreover, analytic) solutions; see a comment at the beginning 
of Section [2] concerning other weaker solutions u > 0. 
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Any rigorous proof of such a blow-up log-log behaviour for (11.241) then becomes rather 
illusive (e.g., more difficult than for the NLSE (jl.6p . and any second-order RD-type prob- 
lem considered before), while formal using of the log-TWs ( 11.1 TP is the only source of a 
proper justification. 

In Section[2]we explain the blow-up log-log phenomenon for the model (11.241) . In Section 
[3j we briefly explain how the same factor occurs in the sixth-order parabolic problem 

(1.25) lit = u 2 (u xxxxxx + u), u = u x = u xx = at x = ±L (L > L ). 

In Sections H] and [5j we show that blow-up is simply self-similar for the better divergent 
models: the fourth- order porous-medium equation with source (the PME-4), 



Without loss of generality, for both PDEs (ll.26|) and (11.27[) . we consider the Cauchy 
problem in R x (0, T) or the IB VP as in (TOD . 

1.8. On extensions to wave and nonlinear dispersion equations. Finally, in Sec- 
tion [HI we show that a similar blow-up log-TW mechanism makes it possible to reconstruct 
log-log factors for some blow-up patterns for other nonlinear PDEs: for the quasilinear 
wave equation (the QWE-4) 

(1.28) u tt = u 2 (-u xxxx + u), 

and for the nonlinear dispersion equation (the NDE-3) 



Actually, the formal matching does not change at all in comparison with the parabolic 
case. Of course, equations (1 1.28ft and (11 . 29[) contain other singularity phenomena such 
as, in view of nonlinear dispersion mechanisms involved (the local speed of propagation 
depends on the value of u itself), formation of various shocks, with rather difficult and 
not fully justified "entropy-like" mathematics (see [23] as a guide for higher-order NDEs 
and [21] for QWEs-4). We do not take these phenomena into account, and just formally 
show that there exist some particular blow-up patterns with log-log property, and do 
not discuss any of their structural stability properties, which can lead to rather obscure 
mathematics. 

1.9. Towards consistency of "topological" matching. As a final comment, the au- 
thor emphasizes that the non-rigorous and rather rough geometric nature of the presented 
results makes no problem for him, since, according to his almost thirty years experience 
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(1.26) 




and for the thin film equation with source (the TFE-4) 
(1.27) u t = ~{u 2 u xxx ) x + u 3 . 



(1.29) 



u t = u 3 (u xxx + u) (note "it 3 " instead of "u 2 "). 



in proving various blow-up results, for some of the above higher-order models^, 



proving the proposed blow-up \/ln In asymptotics 
is not possible in a reasonable finite time. 

In fact, it is not that easy even to explain how difficult those problems are. In Section 
12.51 we present some operator theory comments, and meantime clarify the above point of 
view, partially expressed in (11.301) . in a more formal but consistent manner, as follows: 

(i) It should be noted that the higher-order degenerate singular PDEs such as fll.24p . 
(11.251) . (I1.27p . etc., though having already good and clear applications, have quite obscure 
local existence- uniqueness-entropy-... theory. The problem is not around blow-up (this 
can be done rigorously for some classes of solutions), but even small solutions are extremely 
oscillatory; see similar examples in [T7J and [27J Ch. 3-5]. To say nothing about the 
hyperbolic problem (ll.28p . which cannot have a smooth local solution at all, since in view 
of an obvious nonlinear dispersion mechanism, shock (and rarefaction) waves can appear 
at any suitable point. Note that by no means (I1.28P is a (strictly) hyperbolic system with 
good nowadays entropy theory. Moreover, it can be expected that PDEs such as (11.281) 
principally cannot have a local uniqueness and well-posedness theory. Similar phenomena 
of singularity and non- uniqueness exist for (ll.29p . 

(ii) Overall, in general, in view of a complete absence of local existence-uniqueness and 
entropy theory, studying these higher-order models, we actually deal with huge bundles 
(so-called, flows in dynamical system theory) of solutions, which can be responsible for 
various problems settings, from the Cauchy problem up to infinitely many free-boundary 
problems (FBPs). This is a typical, and unavoidable, feature of modern PDE theory: for 
equations with non-monotone, non-coercive, non- variational... operators of higher-order, 
proving existence, uniqueness, entropy, ... results along the lines great classic theory from 
the twenties century is not only very difficult, but can be impossible in principle. 

(iii) Therefore, dealing with bundles (flows) of blow-up solutions, we propose a rather 
rough asymptotic method, which, as must be admitted, does not specify many par- 
ticular features of solutions involved, but is able to detect some universal property of 
log log-factor, which holds for all such blow-up ones, regardless which problems (Cauchy, 
FBPs, Neumann, Robin, Florin, Stefan,...) are posed. Further refining of the asymp- 
totic method, will require a sharper posing of the Cauchy or FBPs, i.e., matching with 
boundary behaviour, which is completely different (e.g., extremely oscillatory) for various 
problems. Inevitably, one faces difficult local existence/uniqueness/entropy/etc. aspects 
that are unclear for the most of models, which cannot be a subject of a single paper. 

(iv) In a natural sense, the proposed ( "topological" , but very rough, if you wish) match- 
ing well corresponds to common understanding coming from the theory of asymptotic 
series, which are not converging but first terms correctly describe the behaviour of the 
functions relative to a small parameter e — > 0. From planetary motion study in the 
seventeenth century (Newton-Halley's period), it is known that taking many terms of 

8 "The main goal of a mathematician is not proving a theorem, but an effective investigation of the prob- 
lem..." (A.N. Kolmogorov, 1980s; the author apologizes for a non- literal translation from the Russian). 
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(1.30) 



asymptotic expansion gives worse results, while first terms are sufficient to predict the 
behaviour. In a certain sense, a similar asymptotic phenomenon happens in our study: 
trying to improve the expansion, we involve extra, highly oscillatory and singular terms, 
which are deeper connected with the whole solution bundle that we are not unaware of. 

(v) Thus, even in the clear presence of a certain fear of a justified criticism from the 
attentive Readers, the author prefers to keep the style of the asymptotic analysis more 
based on a topological matching, rather than a traditional metric one, where all the terms 
are carefully estimated to convince. Of course, a standard balancing of all the differential 
terms can be done for all the higher-order models, but not more than that: any further 
refining the expansion enters the oscillatory-singular solution bundles and any meaning 
and significance of the blow-up expansion will be lost and become illusive. 

2. HOW LOG-TWS IMPLY THE LOG-LOG BLOW-UP FACTOR 

Thus, we fix the model (ll.24p as the basic one for explaining the occurrence of the log- 
log factor (11. ip in the generic blow-up. Concerning questions of existence and uniqueness 
of solutions of such higher-order equations, we refer to the pioneering paper by Bernis 
and Friedman [7J, where a general approach to constructing such nonnegative solutions 
has been developed. Clearly, their methods of "singular parabolic e-regularizations" for 
constructing non-negative solutions of the Cauchy problem or free-boundary problems 
(depending on parameters) apply to a very wide class of equations including (ll.24p . (I1.25p . 
and many others, linear or nonlinear. On the other hand, for a number of higher-order 
parabolic equations, proper solutions of the Cauchy problem can be oscillatory and of 
changing sign (see [8j [T7J, [27J p. 152] and references therein), but these local interface 
aspects can be neglected at this stage while dealing with peculiar properties of blow-up 
of very large solutions. 

2.1. On nonexistence of blow-up separable-variable solutions. Bearing in mind 
(II. 3p . this is the first crucial step of any blow-up study. Thus, we look for the standard 
self-similar blow-up for (11.241) in the form as in (ll.2ip . 

(2.1) u s (x,t) = ^=0(x), where \ 9 = 9 2 (-9^ + 9). 

Recall that we consider either the Cauchy problem for the equation in (ll.24p . so we 
demand a localized patter satisfying 

(2.2) 9(x) ^ as x -> ±oo, 

or the IBVP as in (11.24=11 . i.e., one needs a proper solution 9(x) such that 

(2.3) 9 = 9' = at x = ±L. 

It is easy to see that such solutions of the ODE in ( 12. II) do not exist. Dividing the ODE 
by 9 2 , multiplying by 9', and integrating by parts yields 

(2.4) | In \9\ = -9"'9 + \ [9" f + \ 9 2 + C, 

where C G IR is a constant of integration. This identity is true in the standard classic 
sense for smooth positive solutions 9(x) > 0, and also remains valid for a proper class of 
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solutions of changing sign having sufficiently regular zeros (this is a question of a suitable 
functional setting, not accented here). We fix this negative result: 

Proposition 2.1. The identity (12. 4p does not allow nontrivial sufficiently smooth, bounded, 
and vanishing solutions 9(x) satisfying ( 12. 21) or f ]2.3j) . 

This negative result can be associated with some king of a strong local monotonicity of 
a nonlinear elliptic operator induced by the equation in (12.11) . 

2.2. On self-similar blow-up for the Neumann or periodic problem. Note that, 
for the Neumann conditions for the IBVP, 

(2.5) 9" = 9"' = at x = ±L, 

the ODE in (12.11) can admit a good proper solution 9(x), and then the blow-up becomes 
purely self-similar as in (12. ip . The same can occur for the periodic problem: 

(2.6) 9{-L) = 0(L), 6'(-L) = 9\L), 9"(-L) = 9"(L), 9"'{-L) = 9"'{L). 

We arrive at the problem of asymptotic structural stability of thast—* T~ , which can be 
solved in a standard manner by considering the rescaled equation for the rescaled solution: 



(2.7) v(x,r) = y/T=tu(x,t), t = -ln(T-t) v T = v 2 {-v xxxx + v) - \ v. 



It is key that (12.71) is a gradient dynamical system with the Lyapunov function 

(2-8) 1 [\ !{v xx f - § / , 2 + \ / In H] = - / %P < 0. 

As in standard blow-up theory, the main difficulty in proving the stabilization as r — > +oo 
in the rescaled PDE (12.71) will be a priori bounds on the orbit, i.e., that 

(2.9) \v(x,t)\<C and sup x \v(x, r)\ -f* as r — > +oo. 

Then, if the set of solution of the ODE (12. ip . (12.51) or (12.61) is discrete, then there exists 
a unique stationary profile 9(x) such that 

(2.10) v(x, t) — > 9(x) as r — > +oo, 

in a suitable metric (e.g., uniformly for strictly positive and hence classic solutions). 

Finally, let us notice that the solvability and multiplicity for the Neumann problem (12.51) 
or a periodic one (12.61) are not easy (though the variational structure suggests existence of 
many periodic solutions), essentially depends on the length parameter L, and can deliver 
some striking surprises. 

2.3. Stationary profiles and a bound on L. Thus, assuming that a standard separable 
blow-up solution does not exist, we then ought to find the stationary one (cf. (11.161) ). 

(2.11) f(x): -f^ + f = on (-L ,L ), f = f = at x = ±L . 
Looking for an even profile, we set (L denotes Lq) 

„ fcos L = — C cosh L, 

(2.12) fix) = cosx + C coshx =^ < . T ~ . , r or tan L = — tanh L. 

[sinL = C sinn L, 

ll 



Matching of log-TW Region-ll with blow-up Inner Region-I 

0.7 1 , , , , 1 , , — 



lnner-1, blow-up 




FIGURE 1. On the idea of a schematic "topological" matching of two Inner Regions. 

The last algebraic equation yields the first positive root L 6 (f,7r), which is taken into 
account in the statement (ll.24p in order to avoid existence of classic stationary solutions 
and hence no blow-up. 

2.4. Formal matching with log-TW behaviour: a "topology" approach. Thus, 
assuming that L > L , similar to (11.141) . we suggest that the blow-up evolution is close 
to the manifold of stationary solutions, i.e., we study the solutions with the behaviour 

(2.13) «CM)~$Mz), 

where f(x) is a stationary solution given by (12.121) normalized by its value at the origin, 

(2.14) /(0) = 1 + C (C > 0). 

The further strategy of matching is explained in Figure [IJ which suggests to match the 
Inner Region-I with the strong and fast blow-up asymptotics (I2.13P with the log-TW Inner 
Region-II situated close to the end point x = L (we consider the x-symmetry geometry). 

We are assuming that in Region-II, the solutions are sufficiently small, so we take into 
account the diffusion-like operator only, 

(2.15) u t = -u 2 u xxxx [x ~ Lq). 
Then, taking the log-TW ffTTTl) yields the following ODE: 

(2.16) \g-Xg' = -g 2 g^. 
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It is now key that the ODE (12.161) admits the following asymptotics: 

(2.17) g(y) = (-y) 2 ^H^)(l + o(l)) as y - -oo. 

Actually, the leading quadratic behaviour ~ y 2 well corresponds to the quadratic "bound- 
ary layer" with the stationary structure 

(2.18) f{x) = C 1 (L -x) 2 (l + o(l)) as x -> L (d > 0). 

The logarithmic correction in (I2.17|) has the same type as in (11.191) . so it is supposed to 
generate the same log-log factor as in (II. ip . 

More precisely, according to (11.171) . we need to match the behaviour in Region-II, 

(2.19) u(x,t) ~ [-x- Aln(T-t)] 2 y/ln(-x - Aln(T-t)), 

with that in blow-up Region-I, where, by the assumption (12.131) . the following holds: 

(2.20) u(x,t) ~ -^Ld(x - L ) 2 for x « L . 

Note that the expansion (I2.19P should be treated as already the extended one from the 
corresponding asymptotic region. A more detailed expansion for smaller values of u will be 
inevitably affected by the singular and oscillatory behaviour of solutions close to x = L 
and x = L, which we would like to avoid at this stage and refer to the discussion in 
Section 11.91 More definitely, the expansion (I2.19P up to some admissible perturbations 
(and partially (12.201) . which is more regular) precisely defines the required class of solutions 
of the IBVP or any of FBP problems. Of course, there are other solutions of the PDE 
under consideration with a completely different blow-up behaviour having nothing to do 
with a double-log factors. 

We omit at this moment the necessary (delicate and questionable) stage in matching, 
where a suitable time- dependence of the log- velocity A = \(t) is assumed, [53, p. 308]. 
Currently, as usual in matching concepts and theory, two manifolds of solutions such as 
(12.191) and (12.201) admit matching if the corresponding pairs of leading multipliers are 
overlapping at some intermediate values of x, possibly, different for both pares. One 
can object that the log-TW expansion in (I2.19P describes moving wave (see the arrow in 
Figured]), unlike the standing one in (12.201) . which, at the first sight, makes matching 
rather suspicious. We should mention in connection with that the log-TW in Region-II 
moves with the logarithmic speed ~ | ln(T — t)\, i.e., very slow in comparison with the 
power rate ~ ^-t °^ ^ e blow-up divergence in Region-I. In other words, relatively, this 
can be classified as an effectively standing wave, which makes the matching correct and 
possible after necessary extra arguments. 

Again carefully looking at the manifolds (I2.19P and (12.201) . we see that the two pairs of 
different multipliers that t -> T- 



(2.21) ( — x - A ln(T - t)) 2 k [x - L c 



and 



— t Vln(-x-Aln(T-t)) k -$L, 



admit natural "structural matching" on some compact subsets (different for both) in x. 
Actually, according to typical concepts of asymptotic analysis, one needs to observe and to 
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get "similar geometric forms" in space and time of matched flows. This can be definitely 
done for the pairs in (12.211) . where the first pair of space distributions are both similar 
and quadratic, while the second pair assumes to be similar for some special functions A(t) 
only. We must admit again that by no means our goal is to prove this matching rigorously 
(whatever "proof means). 
Then setting x ~ L in the second pair leads to the conclusion (II. ip : as t — ■> T _ , 



(2.22) 



^v/ln(-Aln(T-£))~-fL =► A(t) ~ | ln(T - t)|(l + o(l)). 



Note that A even being a slow growing function completely disappears in this matching, 
so even this rough accuracy guarantees the coefficient "1" (and A = X(t) is contained in 
the o(l)-term only, which is not of importance). 

2.5. On some standard metric estimates and operator balance: the origin of 
principal difficulties. Let us return to a more mathematical (rather than geometrical)) 
PDE meaning of these speculations. In fact, this suffices to balance differential terms and 
perturbations in the rescaled equation for v in (12. 7p . where r = — ln(T — t) — > +oo as 
t — > T~ . Actually, Figure [1] is more suitable for the rescaled function v(x,r), which is 
expected to have a slow logarithmic growth as r — > +oo unlike the fast algebraic one for 
u(x,t) as t — > T~. We next need to rescale as in (12.71) the asymptotic flows (I2.19P and 
(1^201) . etc. 

We begin with a standard approach to such a matching by starting with the inner R-I 
rescaled expansion (12.191) . i.e., v(x, r) ~ A(r)f(x), which naturally require introducing 
the new variable w and the corresponding rescaled equation: 

v (x, t) = A(t)w(x, t) and A 2 (r)dr = ds, 



(2.23) 



w s = w 2 (-w xxxx + w) + irkr-K w - w - 



2A 2 (t) w A 3 (t) 

Thus, here A{t) is unknown, but (II. ip suggests that 

(2.24) A(t) ~ \f\ nr and s~rlnr (r ~ ^) as r — > +oo (expectation). 

Indeed, the perturbed equation ( 12.231) with the leading (by ( II. ip ) perturbation ~ O ( ) 
describes the stabilization to the stationary profile: 

(2.25) w(x,r) — > f(x) as r — >• +oo, where f(x) solves (12. lip . 

on compact subsets in [0, L ). Note that, for the desired result (I2.24p . the following holds: 
(2-26) ^^^L-^LK^oo)), since ~ dr, 

i.e., the perturbation is not integrable in L 1 that causes extra difficulties (though is a 
typical feature of many blow-up and extinction problems, [28]). However, this stage of 
the analysis is rather standard, since the unperturbed equation 

(2.27) w s = w 2 (-w xxxx + w) 

is a gradient system (cf. (12.81) ). though fighting the perturbation terms in (I2.23P can be 
rather technical (but hopefully not of a principle issue). In standard blow-up theory such 
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problems are tackled by theorems on stability of u;-limit sets with respect to arbitrary 
small perturbations of the dynamical system, provided that the unperturbed problem is 
uniformly Lyapunov stable; see examples and a large number of references in [28] . 

Anyway, this is not the case, where the main difficulty occur. Namely, we need to 
match the behaviour (I2.23p . (I2.25|) with that corresponding to a proper one induced by 
the boundary conditions at x = L (or an FBP setting, which is hard to distinguish from). 
Then we ought to perform the linearization: 

(2.28) w(x, s) = f(x) + Y(x, s) =s- Y s = CY + 2fY(-Y^ + Y) + ^ f + ... , 

where we omit all asymptotically smaller linear and nonlinear terms, which however can 
be important for matching. As usual, equation in (I2.28P clearly shows that full spectral 
theory of the linearized operator 

(2.29) C = f(x)(-D 4 x + I) in x G [0, L Q ) 

is necessary for extension of the solution for x > L . Fortunately, it is symmetric in the 
weighted space L 2 , with the weight p = f~ 2 (x) > 0, which is singular at the end point 
x = L . However, checking the crucial end point x = L Q and using the asymptotics of the 
coefficient ( 12.181) . yields Euler's equation in the spectral problem: 

(2.30) C 2 (L - x) 4 (-^ (4) + V>) + - = M>- 

Hence, linearly independent solutions polynomial with the characteristic equation: 

(2.31) ip(x) ~ (L - x) m Cfm(m - l)(m — 2)(m - 3) + A = 0. 

This shows that the necessary spectral theory, though covered by classic theory of singular 
ordinary differential operators (see, e.g., Naimark |38j). can be rather tricky. In particular, 
deficiency indices of the operator (12.291) are not easy to detect in L 2 More precisely, 
depending on the boundary matching, this can require self-adjoint extensions in both 
the cases of discrete and continuous spectra. Infinitely oscillatory behaviour for complex 
values ofm = a + i&GCin (12.31 j) (A remains real) occurs close to Lq , where 

(2.32) Y(x, s) ~ e Xs (L Q - x) a sin(Mn(L - x)) as s — ► +oo: 

cf. the discussion of oscillatory issues in Section 11.91 Thus, using proper spectra theory 
of (12.291) with eigenfunction expansions of w(x,s) of (I2.28P over a discrete (better) or 
continuous (most plausible for many problem settings) is the first key concern here. 

Continuing the spectral aspects of the expansion, we next assume that the spectrum of 
C in (12.291) is continuous in a suitable setting and deficiency indices, so we may impose 
the following behaviour of the solution Y of (12. 28ft : 

(2.33) Y(x, s) = fi(s)$(x) + ... as s -> +oo. 
Substituting this into (I2.28P and naturally assuming that 

(2.34) = = =► \»'(s)\ « \fx(s)\ for s»l 

(this is true for the expectation (I2.24p ). we arrive at the following first rough balance: 

(2.35) // s $ = /iC$ + ^/ + ... => M«) = 2A% and C$ + / = 0. 
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Again, a suitable solvability of the ^-equation in f!2.35j) requires posing proper boundary- 
conditions at the singular end-point x = L and satisfactory spectral theory of C. 
Meantime, we solve the ODE for $ in (12.351) close to the singular end-point x = Lq: 

(2.36) -*« + $ = = -a^^ as x - L . 
This yields the following non-regular asymptotic expansion: 

(2.37) $(x) = z 2 \nz + A l + A 2 z + A 3 z 2 + A 4 z 3 + ... as z = L -x-^0 + . 

Here we see the first clear appearance of the ln(L — x) in the framework of this standard 
matching and operator analysis: 

(2.38) $(x) ~ (L Q - x) 2 ln(L - x) + A x + A 2 {L - x) + ... as ar -> L , 

which is indeed a reliable remnant of the asymptotic expansion (I2.17P for the log-TW 
profile g(y) satisfying (12.161) . The integration constants ^2,3,4 G K in (12.371) should be 
obtained by further matching (some of them can be arbitrary and depend on initial data). 

Next, in the R-II, where w is uniformly small and hence w 3 can be neglected, the 
problem (I2.23f) . with the behaviour (I2.33H at the end-point x = Lq, takes the form 

w s = -w 2 w xxxx + fi(s)w + ... in (L ,L)xR + , 
w = Axfi(s), w x = —A2fJ,(s), x = L ; w = w x = 0, x = L. 



(2.39) 



The precise behaviour of w and the derivatives w x and w xx at x = L is defined by (I2.33|) . 
Therefore, we add to (I2.39P the following overdetermining condition 

(2.40) w xx (x, s) = — ln(x — Lq)(1 + o(l))/i(s) as x — > Lq , 

which is expected to detect the unknown function A(r). Theoretically, once these have 
been achieved, extending the solution up to the original boundary point x = L will define 
the necessary ODE-type equation for the unknown function A(r) such as, very roughly, 

(2.41) A' ~ e~ A2 A(t) ~ v^nr for r > 1. 

The principle key difficulties start now, when we arrive at the problem (12.391) . The 
point is that when w(s,x) approaches 0, the problem becomes extremely sensitive and 
demand solving a number of linear and nonlinear eigenvalue problems, where we face in 
a full scale those unresolved local existence-uniqueness-entropy-etc. problems mentioned 
in Section 11.91 Namely, even the unperturbed problem 

(2.42) w s = —w 2 w xxxx on small solutions \w\ 1 

is not well-understood at all. Note that the second-order counterpart of (ll.lOp . 

(2.43) w s = w 2 w xx , 

is also not easy, but by a contact Backlund symmetry reduces to the heat equation: 

(2.44) v s = v vv , where w(x,s) = -r-r — r and <^ s 

y ' s yy, \ , ) <f> x {x,s) ]x = V{y,S); 
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see [53J P- 78] and references therein. Actually, such a homology (I2.43[) i— > (I2.44p eventually 
creates a relation (12.411) . where the exponential term e~ A is naturally associated with the 
Gaussian fundamental solution of the operator D s — Dy in (12.441) . Obviously, no such 
nonlocal symmetry is available for the much more difficult PDE (I2.42p . 

Thus, the technicalities at this last stage are extremal and which deserves extra special 
difficult analysis that we cannot be performed in the present more formal paper (will 
require at least an extra dozen of sheets of hard calculations), so that, mathematically: 

Main difficulty: problems (12.391) and (I2.42p are not doable at the moment 

(an open problem for the attentive Reader). Let us note that spectral theory of generalized 
fourth-order non-symmetric Hermite operators (an operator pair) [15] can be key: 



As we have mentioned, exactly in this R-II, a delicate and specially deformed structure 
of log-TWs are crucial for creating the whole pattern. Thus, the geometric approach of 
formal matching makes it possible to avoid some (too) detailed refinement of the solution 
expansion, and establishes that the double-log factor can occur for a kind of an "envelope" 
of a wide class of blow-up solutions induced by various functional settings. However, 
the behaviour is expected to be so difficult that we will not be surprised (e.g.,) if the 
resulting more correct calculation provides us with a different double-log factor A(t) ~ 

[in | ln(T — t)|] 3 ^ 4 instead of the standard one (11.11) . though the topology matching via 
(12.211) . (I2.22p indicates that this is not the case. 

Finally, let us point out again that a full (formal) asymptotic expansion technique 
also assumes the next matching of the log-TW Region-II with the Boundary one close 
to x = L, where another (and hopefully and typically much weaker than at x = L ) 
boundary layer of uniformly bounded solutions can occur. As usual, this is expected to be 
easier (and less principal), but can be also tricky, as everything concerning such nonlinear 
degenerate higher-order operators. We recall that a proper and mathematically justified 
matching of those manifolds of solutions via rigorous compactness-like arguments seems 
not achievable for the higher-order parabolic equations under consideration. Note again 
that the analysis in [1] for (11.101) was essentially and inevitably based, among others, on 
Maximum and Comparison Principle arguments, which are non-existent for higher-order 
parabolic flows. Moreover, as we have pointed out in greater detail in Section 11.91 the 
CP, IBVP, or any FBP settings for (ll.24p have not been well developed still, so we have 
to be very careful in deriving asymptotic blow-up properties of some possibly nonexistent 
and/or non- unique solutions. 

3. Briefly on the sixth-order model: the log-log universality 



(2.45) 




!, k>0}. 



For the equation ffl~25j) . the ODE in (TJH]) is 



(3.1) 
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so that nonexistence of a sufficiently smooth solution follows from the identity (cf. (12.41) ) 

(3.2) \ In \Q\ = 0^6' - 0^6" + \ {Q"f + \6 2 + C. 
The ODE for the log-TWs in Inner Region-II now takes the form 

(3.3) \g = #V 6) + - #0?) ~ i73 (-^) 3 \/M Z ^) as r)^-oo. 

Therefore, the same argument of matching with the stationary structure in blow-up Inner 
Region-I governed by 

(3.4) f(x) : + f = in (-L , L ), / = /' = /" = at x = ±L 

(hence with the cubic behaviour f(x) ~ Ci(L — x) 3 for x ~ L that is well matched with 
the term (— rj) 3 in ( 13. 31) ) yields as in (12.22ft the same log-log factor A(t). 

4. NO LOG-LOG FOR THE PME-4 WITH SOURCE 

4.1. Countable basic family of patterns by variational approach. Consider the 
Cauchy problem for (11.261) in R x (0, T). Following the question (ll.3p . we then claim that 
there exist weak compactly supported similarity solutions in the separable variables: 

(4.1) Us (x,t) = ^= t 6(x), where \ 6 = -(9 3 )^ + 9 3 in R. 
Changing for convenience the function, 

(4.2) 0(x) = j s Fi(x), 

yields a more standard semilinear ODE with non-Lipschitz nonlinearity, 

(4.3) F (4) - F + = in R, 

which admit three obvious constant equilibria and ±1. Fortunately, this problem is vari- 
ational, and the following result is proved by a combination of Lusternik-Schnirel'man 
category (genus) theory (L-S theory, for short) and Pohozaev's fibering approach of cal- 
culus of variations (26]. 

Proposition 4.1. The problem (14.31) admits at least a countable set of nontrivial com- 
pactly supported solutions. 

We present a few comments concerning this principal result. We look for critical points 
of the corresponding ^-functional: 

(4.4) E{F) = -\J{F»Y + \$F*-\j\F\t 

In general, we have to look for critical points in W$(M. N ) R L 2 (R N ) n L^ 3 (R N ). Bearing 
in mind compactly supported solutions, we choose a sufficiently large R > and consider 
the variational problem for ( 14. 4ft in W$ q (Br), Br = (-R, R), where we assume Dirichlet 
boundary conditions. It is next proved that any solutions satisfying F(y) — > as \y\ — > oo 
is compactly supported in R. 

Thus, the functional (14.41) is C 1 and is uniformly differentiable and weakly continuous, 
so we can apply classic Lusternik-Schnirel'man (L-S) theory of calculus of variations [3U 
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§ 57] in the form of the fibering method [501 EE]. Namely, following L-S theory and the 
fibering approach [5T], the number of critical points of the functional (14.41) depends on 
the category (or genus) of the functional subset on which fibering is taking place. Critical 
points of E(F) are obtained by spherical fibering 

(4.5) F = r{v)v (r>0), 

where r(v) is a scalar functional, and v belongs to a subset in W$q(Br) given as follows: 

(4.6) H = {v G Wl (B R ) : H (v) = - j\v"f + J = l}. 
The new functional 

(4.7) H(r,v) = \r 2 - § A J 
has the absolute minimum point, where 

(4.8) H' r = r-A J\v\l =0 r Q (v) = (J \v\*)* . 
We then obtain the following functional: 

(4.9) H(v) = H(r (v),v) = -±r » = -\ {J \v\lf. 

The critical points of the functional (14. 9 ft on the set (14.61) coincide with those for 

(4.10) H(v) = J \v\l, 

so we arrive at even, non-negative, convex, and uniformly different iable functional, to 
which L-S theory applies, j3H § 57]; see also [HI p. 353]. Following [5T], searching for 
critical points of H on the set Ho one needs to estimate the category p of the set TLq. The 
details on this notation and basic results can be found in Berger [6], p. 378]. 

It follows that, by this variational construction, F is an eigenfunction satisfying — F^ + 
F — fiFa = 0, where /i > is Lagrange's multiplier. Then scaling F i— > /x 3//2 F yields the 
original equation in (14.31) . 

For further discussion of geometric shapes of patterns, it is convenient to recall that 
utilizing Berger's version [61 p. 368] of this minimax analysis of L-S category theory [3U 
p. 387], the critical values {ck} and the corresponding critical points {vk} are given by 

(4.11) c k = 'mf^ Mk sup v6JF H(v), 

where T C 7io are closed sets, and M.^ denotes the set of all subsets of the form BS^ 1 C 
TCo, where S k ~ l is a suitable sufficiently smooth (k — l)-dimensional manifold (say, sphere) 
in 7i and B is an odd continuous map. Then each member of Aik is of genus at least k 
(available in Hq). It is also important to remind that the definition of genus [2U p. 385] 
assumes that p{J-) = 1, if no component of T U J 7 *, where T* = {v : v* = —v G J 7 }, 
is the reflection of T relative to 0, contains a pair of antipodal points v and v* = —v. 
Furthermore, p(^ r ) = n if each compact subset of T can be covered by, minimum, n sets 
of genus one. 
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According to (14.111) . c\ < c 2 < ... < q , where l = lo(R) is the category of TCq (see an 
estimate below) such that 

(4.12) l (R) -> +00 as R^oo. 

Roughly speaking, since the dimension of the sets JF involved in the construction of M. k 
increases with k, this guarantees that the critical points delivering critical values (14.111) 
are all different. 

It follows from (14.61) that the category l = p(7i ) of the set Ho is equal to the number 
(with multiplicities) of the eigenvalues X k > —1, 

(4.13) l = p (H ) = ${\ k > -1} 
of the linear bi-harmonic operator — D A < 0, 

(4-14) - V>£ 4) = \ k i>k, ipk e W! >0 (B R ); 

see J6j p. 368]. Since the dependence of the spectrum on R is, obviously, 

(4.15) X k (R) = R- 4 X k (l), k = 0, 1,2, ... , 

we have that the category p(H.o) can be arbitrarily large for R ^> 1, and (I4.12p holds. 

4.2. On total variety of patterns: numerical evidence. Actually, the total variety of 
possible solutions is not exhausted by the categories of L-S theory, which we will explain 
below by presenting clear numerical evidence. 

Figure [2] demonstrates the first, and actually, the ground state profile Fq(x). It has a 
typical shape of a ground state, as a critical point of the functional delivering the absolute 
extremum. Note that unlike the classic second-order case of the ground state for [12J (this 
problem is key for the critical NLSE) 

(4.16) Am-« + m 3 = in R N , 

which is strictly positive (with exponential decay at infinity) and is unique up to transla- 
tions, the ground state F for (14. 3p is oscillatory and of changing sign near finite interfaces; 
this is seen from Figure [2j The oscillatory structure of solutions is rather involved and is 
described in [26] ; see also similar details in [T7] . 

In Figure[3j we show the second dipole-like profiles, where F±(x) denoted by the boldface 
line is the basic one corresponding to the L-S critical value. In addition, there exists a 
countable (this is not proved still) family of dipoles {F^}, which differ by the structure 
of the internal zero close to the origin (in general, k stands for the number of transversal 
zeros there, but this is not enough to uniquely identifying the pattern). By the dotted 
line, we denote another profile from the next family {i^ }i see below. 

Figure H] continues explaining further basic L-S patterns, where we show F , Fi, F 2 , and 
F3 (the dotted line). It is clearly seen that each F k has precisely k "dominant" transversal 
zeros inside the support, which well corresponds to Sturm-like principle (not applicable 
here in the rigorous sense, since all the solutions are oscillatory and have infinitely many 
sign changes near interfaces). 
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The firs (ground state) profile F (x) 




FIGURE 2. The first, ground state solution of the ODE (|43 



Dipole— like profiles F* k '(x 
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FIGURE 3. Some of the dipole-like solutions F^\x), k = 0, 1,2 of the ODE (TO) . 



In Figure [5j non L-S profiles from the family {F + 2k} are presented. Figure [6] shows 
some profiles from the family {F +2i 2fc,+2(^)}, which are also not expected to correspond 
to any L-S category/critical value. 

First profiles from the family {i 7 +2,2,...,2,+2(a ; )} (also non L-S) are shown in Figure [71 
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First four basic (L-S) profiles 
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FIGURE 4. Four basic L-S patterns of the ODE <|43j) 
First two profiles from family {F } 
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FIGURE 5. Solutions F +4 (x) and F +6 (s) of the ODE (Oil . 



Finally, following [26] and using the above rather simple families of patterns, we claim 
that a pattern (possibly, a class of patterns) with an arbitrary multiindex of any length 



(4.17) 



a = {±cri, 0- 2 , ±a 3 , 04, ±07} 

22 




-20 -15 -10 -5 5 10 15 20 

X 

FIGURE 6. Solutions {F +2 , 2 k,+2(x)} of the ODE (jOj) . 



F(x)i- 
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X 

FIGURE 7. Solutions {i ? + 2 ,2,...,2,+2(x)} of the ODE (|C5|h 

can be constructed. Here, as above, each itr^ stands for the total number of succes- 
sive intersections with the current equilibrium ±1, while counts that with the trivial 
equilibrium 0. E.g., in Figure [H we show a single complicated profile F a [x) with the index 



(4.18) 



a = {+4, 1, -4, 1, +2, 2, +4, 1, -4, 1, +8}. 
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Complicated pattern F o with o={+4,1 -4,1 ,+2,2+4,1,-4,1 ,+8} 




-50 50 

X 



Figure 8. Solution 1—4 1+2,2, +4,1,— 4,i,+8 

0) of the ODE (TOl) . 

Actually, the multiindex (14. 1 7f) can be rather arbitrary (with some natural restrictions 
on even local intersection numbers) and then takes finite parts of any admissible non- 
periodic fraction. Overall, this means chaotic features of the whole family of solutions 
{F a }. These chaotic types of behaviour are known for other fourth-order ODEs with 
coercive smooth operators, [45, p. 198]. 

5. NO LOG-LOG FOR THE TFE-4 WITH SOURCE 

For the TFE-4 (ll.27p . the similarity substitution as in (14.11) yields another ODE: 

(5.1) \e = -(e 2 0'"y + e 3 in r. 

To get equilibria ±1, we perform the change: 

(5.2) 6{x) = X F(x) (F 2 F"')' = F 3 — F in R. 

Unlike the one in (14.11) and (14.31) . the ODE in (15. 2p does not admit a variational formula- 
tion. Nevertheless, there exist rather standard shooting arguments for detecting necessary 
similarity blow-up profiles; see references in [T7J §3]. In particular, the following asymp- 
totic behaviour is known close to finite interfacesj 

(5.3) F(x) = \ (xq — x) 2 + [a/| ln(x — x)\ + C + ...] as x — > x , 

where C G R is an arbitrary constant. Overall, the asymptotic bundle ( 15. 3D comprises 
two parameters {xq,C}, which are expected to be sufficient to shoot also two conditions 



'These are the maximal regularity solutions admitted by the ODE (|5.2[k see [T7] . 
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First even pattern F (x): x =2.83... 



F(x) 

0.5 









x =2.83... 

v/ 
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X 



FIGURE 9. The first basic nonnegative solution Fq(x) of the ODE in (|5.2I) , 
at the origin: 

(5.4) F(0) = F"(0) = (symmetry) or F(0) = F"(0) = (anti-symmetry). 

Figure [9] shows the first nonnegative even profile Fq(x) for x > with the symmetry 
conditions in ( 15.4ft and the interface xn = 2.83... . Hence, the answer to ( 11 .3D is log-log jB. 

6. QUASILINEAR WAVE AND NONLINEAR DISPERSION EQUATIONS 

6.1. The QWE— 4. For (11.281) . the similarity separable solution (12.11) is slightly different, 
but the ODE remains the same: 

(6.1) u s {x,t) = ^r t 9(x), where 29 = 9 2 (-9^ + 9). 
Hence, Proposition 12.11 holds, so we need the log-TWs for matching: 

(6.2) u log (x,t) = ^- t g( V ), V = x + X\n(T-t) => 2g - 3Xg' + X 2 g" = g 2 (-g^ + g). 
However, the necessary asymptotics is guaranteed by keeping just two terms: 

(6.3) 2 g = -g 2 g^ + ... : 

so that the only correction to ( 12. 17ft is the multiplier a/2. The rest of the analysis remains 
the same as at the end of Section [2j 
For the divergent wave model 

(6.4) u tt = -{u 3 ) xxxx + u 3 

the log- log factor is nonexistent, since it admits separable compactly supported solutions 

(6.5) u s {x,t) = ^r t 9{x) 29 = -(9 3 )^ + 9 3 , 

25 



which on scaling reduces to that in (14.11) . 

6.2. The NDE— 4: non-divergent model. For (11.291) . the similarity separable solution 
(12. ip is slightly different, but the ODE remains the same: 

(6.6) u s (x,t) = (T-t)-*9(x), where § 9 = 6 S (6"' + 9). 

Instead of (12 .4p . the nonexistence conclusion is governed by the following "local mono- 
tonicity" identity: 

(6.7) -|| = 9'9" - J(9") 2 + \e\ 

so that sufficiently smooth solutions with 9 = at ±L (and 9 e H 2 ) do not exist. 
As the analogy to (I2.12p . consider the stationary equation: 

(6.8) f(x): f"' + f = on (-L ,L ), f(-L ) = f(-L ) = 0, f(L o ) = 0, 

where the zero at x = L is assumed to be transversal with the behaviour (cf. the 
smoother one (I2.18P ) 

(6.9) f(x) = d(x - L )(l + o(l)) as x -> L {C x > 0). 

According to the problem (16.81) . to avoid extra difficulties with rather obscure conse- 
quences, we consider (I1.29P in (— L , L) x (0, T), with L > L (possibly, L ^> L to avoid 
existence of other stationary profiles) and the same boundary conditions at the left-hand 
end point, 

(6.10) u = u x = at x = —L and u = at x = L. 

Therefore, the log-log perturbation will penetrate into Inner Region-I from the right-hand 
singular point x = L , where the stationary solution (16. 8p vanishes thus creating an 
internal singular layer to be resolved by using slowly moving log-TWs. 
Thus, the log-TW ansatz yields the ODE 

(6.11) U log {x,t) = (T-tyh(v), V = X + XHT~t) ^lg-\gi = g 3 (g'» + g). 

Close to the necessary point x = L , we keep as usual two terms that yields the desired 
behaviour: 

(6.12) \g = g 3 g"' + ... =^ g(r}) ~ (-77) [In (-??)] 3 as 77^-00. 

Similar to (I2.21I) - (I2.22I) . matching the linear structure (—77) in (16.121) with the linear one 
~ (L — x) in (16.91) yields 

(6.13) A{t) ~ [ln|ln(T-f)|p(l + o(l)) as t^T~. 

Observe changing y(-) into the cubic root due to cubic nonlinearity u 3 in (ll.29p instead 
of u 2 in other models. The choice of the quartic equation (I1.29P was in fact generated by 
the possibility of such a transversal matching of "linear" structures. Some aspects of such 
a matching are still unclear and deserve further study. 
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Typical blow-up similarity profiles F(x) for the NDE 




FIGURE 10. Blow-up patterns F{x) satisfying (|6.16p with the behaviour (|6.17p . 
6.3. No double-log in a divergent NDE. For the corresponding divergent NDE, 

(6.14) u t = {u A ) xxx + u 4 , 
answering (11.31) . one can construct separable blow-up solutions 

(6.15) u s {x,t) = {T -t)-U{x) | 9 = (6> 4 )"' + e\ 
Looking for nonnegative solutions 9 > 0, on rescaling, yields 



(6.16) 



9 = 3"F 



F"' + F-F*=0 in (-L ,L oj 



Unlike the ODE (I4.3p . this one admits solutions with a single interface (free-boundary) 
point x = — L with the behaviour 



(6.17) 



F(x) = (24) ~s (Lq + x)Ul + o(l)) as x 



-Li 



In other words, we can set F(x) = for all x < —Lq. Therefore, moving such a profile 
F(x), one can always satisfy the boundary condition at x = L in ( 16. 8ft for any L > 0. 
Typical solutions of the ODE (I6.16P are shown in Figure [10l where positive humps can 
serve as blow-up patterns. Note that all such solutions are not oscillatory near finite 
interfaces, unlike those for the TFEs [17], where another third-order oscillatory ODE 
occurs. Thus, the similarity law (I6.15P describes blow-up in the divergent model (I6.14p . 



7. Final conclusion: blow-up log-log is universal in PDE theory 
The above study allows us to fix the following conclusion: rather surprisingly, 

(7.1) 



the factor yln | ln(T — t)\ has a clear universality in blow-up 
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for different classes of higher-order (non-divergent) nonlinear evolution PDEs. The In ln- 
factor also occurs in boundary regularity (Petrovskii-type) analysis for equations [22] 

1H U X xxxi U'xxxi V>t )xrau etc. 

It would be important, on the basis on the above discussion of various linear and nonlinear 
PDEs and by adding new asymptotic phenomena of necessity, to explain, in a more 
unified way, and to derive, formally or more justified, such a common matched asymptotic 
Vlog log-criterion. 
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